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To clarify the origin of a charge-density wave (CDW) phase in ir-TiSe 2 , we study the ground 
state property of a half-filled two-band Hubbard model in a triangular lattice including electron- 
phonon interaction. By using the variational Monte Carlo method, the electronic and lattice degrees 
of freedom are both treated quantum mechanically on an equal footing beyond the mean-field 
approximation. We find that the cooperation between Coulomb interaction and electron-phonon 
interaction is essential to induce the CDW phase. We show that the “pure” exciton condensation 
without lattice distortion is difficult to realize under the poor nesting condition of the underlying 
Fermi surface. Furthermore, by systematically calculating the momentum resolved hybridization 
between the two bands, we examine the character of electron-hole pairing from the viewpoint of 
BCS-BEC crossover within the CDW phase and find that the strong-coupling BEC-like pairing 
dominates. We therefore propose that the CDW phase observed in ir-TiSe 2 originates from a 
BEC-like electron-hole pairing. 

PACS numbers: 71.10.-w, 71.45.Lr, 71.35.Lk, 71.27.-|-a 


I. INTRODUCTION 

Charge-density wave (CDW) is widely observed in low¬ 
dimensional solids and has been extensively studied both 
experimentally and theoretically [I|, Q . Transition metal 
dichalcogenides MX 2 (M=transition metal, X=S, Se, Te) 
are one of the typical CDW materials with a layered 
triangular lattice structure. They show various CDW 
patterns, depending on the combination of M and X 
and quite often superconductivity (SC) is observed 
next to the CDW phase by applying pressure dop¬ 
ing ii, or intercalation [T3|- However, the origin of 
CDW and SC has not been fully understood and it has 
been still under debate in spite of the long and extensive 
studies so far. 

Recently, lT-TiSe 2 , one of the old transition metal 
dichalcogenides, has again attracted much interests in 
the context of exciton condensation. This material is a 
semimetal or a semiconductor in room temperature [n^ 
[T^ and shows a commensurate CDW transition with a 
2x2x2 superstructure below Tc ^ 200K. The Fermi sur¬ 
face (FS) partially remains even in the CDW phase due 
to the imperfect opening of the energy gap. The origin of 
the CDW phase is still controversial and the usual nesting 
mechanism seems unlikely due to its poor FS nesting M- 
Alternatively, exciton condensation has been propos ed as 
a possible mechanism for the CDW phase |15l - [r^ . In¬ 
deed, the large spectral weight transfer between Ti 3d 
and Se 4p bands and the flat energy spectrum just below 
the Fermi energy have been observed [l^, strongly sug- 
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gesting the possibility of exciton condensation. On the 
other hand, another possible mechanism for the CDW 
transition is a band Jahn-Teller effect which results from 
electron-phonon interaction [Tsl - I^ . A large lattice dis¬ 
tortion of several percent observed below Tc 0 indicates 
the strong coupling between electronic and lattice de¬ 
grees of freedom [22l| . Very recently, it is proposed that 
both mechanisms work cooperatively for the CDW tran¬ 
sition (23 I - It^ . 

Exciton condensation is a quantum state expected in 
low carrier density systems such as a semimetal or a 
semiconductor and has been extensively studied since 
1960s (2§ - [^ . The exciton is a bound pair of an electron 
and a hole in different bands mediated by the interband 
Coulomb interaction. When the binding energy of an 
exciton exceeds the band gap, the system has an insta¬ 
bility toward condensation of excitons, namely, a spon¬ 
taneous hybridization between different bands. Since the 
repulsive Coulomb interaction is attractive between an 
electron and a hole, the exciton condensation is expected 
to occur in principle without considering any additional 
“glue” of the electron-hole pair. Although extensive ef¬ 
forts have been devoted for half a century and several 
candidates have been proposed [B IB - lssj l. the gener¬ 
ally accepted materials for the exciton condensation are 
still absent. Therefore, any conclusive evidence for the 
exciton condensation in real materials is highly desired 
for further progress and lT-TiSe 2 would be one of the 
promising examples. 

In the pioneering studies for exciton condensation [B . 
[B| . an isotropic band dispersion with perfect FS nest¬ 
ing is assumed, for which the excitonic instability is al¬ 
ways present in a semimetallic case. The extension to 
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anisotropic band dispersions shows that the degree of FS 
nesting greatly affects the instability of exciton condensa¬ 
tion [3^. Moreover, the mean-field approximation gener¬ 
ally overestimates the instability toward ordered states, 
including the exciton condensation. Therefore, for dis¬ 
cussion of the exciton condensation in real materials, the 
realistic band dispersion and the appropriate method be¬ 
yond the mean-field approximation are both required. 

In this paper, a half-filled two-band Hubbard model 
with electron-phonon interaction in a triangular lattice 
is studied to understand the origin of the CDW phase 
in lT-TiSe 2 . The ground state properties are calculated 
using the variational Monte Carlo (VMC) method for 
multiorbital systems [s^- We find that the cooperation 
between Coulomb interaction and electron-phonon inter¬ 
action is essential to induce the CDW phase. The CDW 
phase is observed between the normal metal and band 
insulator phases with intermediate interband Coulomb 
interaction. We show that the “pure” exciton condensa¬ 
tion without lattice distortion is difficult to realize un¬ 
der the poor FS nesting condition in a triangular lattice. 
Furthermore, we systematically calculate the momentum 
resolved hybridization between the two bands to show 
that the strong-coupling BEC-like pairing dominates in 
the CDW phase. Our results therefore suggest that the 
CDW phase observed in lT-TiSe 2 originates from the 
strong-coupling BEC-like electron-hole pairing. 

The rest of this paper is organized as follows. In Sec.lBl 
a two-band Hubbard model in a two-dimensional trian¬ 
gular lattice is introduced as a low energy effective model 
for lT-TiSe 2 . The detailed explanation of the VMC 
method and the variational wave functions are also given 
in Sec. m The numerical results are then provided in 
Sec. IBII Finally, the implication of our results for IT- 
TiSe 2 is discussed in Sec. IIVI followed by the summary 
in Sec. [V] 


II. MODEL AND METHOD 


We consider a two-band Hubbard model in a two- 
dimensional triangular lattice defined as 


^ ~ + '^^kfkafkcr 

k,(T k,<7 

i i i 

^ k^q,(T 

q ^ / 

where c\.^ ifta) creates an electron in c (/) band with 
momentum k and spin (7(=t,4,). The band dispersions 



FIG. 1. (color online) (a) Fermi surfaces (blue circle and 
red ellipses) and (b) energy dispersions in the noninteract¬ 
ing limit for the two-band Hubbard model with electron 
density n = 2. A set of tight-binding parameters used is 
{tf, tc, t'c, He) = (1.0, 0.2, 0.15, 6.0) t with t = tf as an energy 
unit, (c) Schematic real-space figure of a 2x2 superstructure 
(dotted lines) with each supercell containing 4 sites (A, B, 
C, and D). The corresponding ordering wave vectors qi, q2, 
and qs are indicated in (a). Dashed lines in (a) represent the 
folded Brillouin zone due to the formation of the 2x2 super¬ 
structure in (c). (d) Schematic real-space figure of the lattice 
distortion considered in the c-f plane along one particular di¬ 
rection, e.g., AB direction indicated in (c). m represents the 
displacement of c atom at site i from its original position. 


and si. are given as 


= 2tc ^cos kx + 2 cos ^kx cos J 

/ 0 ^ \ 

-I- 2t'^ I cos V^ky + 2 cos -kx cos 1 + Me (2) 


and 


el = 2tf (coskx+2cos^kxCOS^ky\ , (3) 


respectively. We introduce the next-nearest-neighbor 
hopping to locate the bottom of c band at M points. Uc 
(D/) is an on-site intraband Coulomb interaction within c 
(/) band and U' is an on-site interband Coulomb interac¬ 
tion between c and / bands, nf^ is a number operator of 
a (= c, /) electron at site i with spin a and 
q(fc, q) is an electron-phonon coupling constant and is 
a creation operator of phonon with momentum q and fre¬ 
quency uj{q). In this model, the lattice distortion changes 
the c-f bond length as shown in Fig.jljd), which couples 
to the c-f hybridization modulated with wave vector q 
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through g{k,q). The total number of sites is indicated 
by N. 

The noninteracting tight-binding parameters, 

{tf, 4, t'„ He) = (1.0, 0.2, 0.15,6.0) t, (4) 

are set to mimic the electronic structure of lT-TiSe 2 with 
a hole pocket at the T point and electron pockets at the M 
points as shown in Figs.[TJa) and[TKb). The ordering wave 
vectors connecting F and M points are denoted as qi, q 2 , 
and qa. Although the ordering wave vectors observed in 
lT-TiSe 2 connect F and L points with a finite compo¬ 
nent , here we consider a pure two-dimensional model 
for simplicity and the limitation of the model is discussed 
later. 

The effect of Coulomb interaction and electron-phonon 
interaction is treated on an equal footing using a VMC 
method. We consider the trial wave function as follows: 

|vk)=Pe-ph|'I'ph)|4/e). (5) 

|4>e) = Pq Pjc 1 ^) is an electron wave function con¬ 
sisting of three parts. |d>) is a Slater determinant con¬ 
structed by diagonalizing the one-body part of Hamilto¬ 
nian H including the variational tight-binding parame¬ 
ters {tf = 1, tc,t'^,Hc) and the off-diagonal element V 
which induces the c-f hybridization. Here, we assume 
V = Vb exp[—and Vq and A are both 
variational parameters. Vq is an amplitude of the c-f 
hybridization and A controls the internal extent of ex- 
citon in k space. We have found that the variational 
energy is improved by introducing A and the behavior of 
A is related to the BCS-BEC crossover of exciton con¬ 
densation [1^. is a Gutzwiller factor extended for 

two-band systems [11,14^. In Pq \ possible 16 patterns 
of charge and spin configuration at each site jF), i.e., 
|0) = |0 0), |1) = |0 t), • • •, |15) = Iti U), are differently 
weighted and their weight {gr} are optimized as varia¬ 
tional parameters. = exp[- Hap is 

a charge Jastrow factor which controls long-range charge 
correlations. Here, = v°‘^{\ri — rj\) is assumed and 
Vi is the position of site i. 

The trial wave function for phonon is assumed to be 
a Gaussian in the normal coordinate {Qq} representa¬ 
tion [4ll.[43|. 

= ({QJl^-P*') = exp 

where Qq = /y/N is Fourier transform of real 

space lattice distortion {rti} at site i. Since the ordering 
wave vectors qi, q 2 , and qa are exactly half of the recip¬ 
rocal lattice vectors of the normal phase [see Fig. [Il[a)], 
the corresponding normal coordinate Qq^ (i = 1,2,3) 
are real numbers. Therefore, we can take the trial wave 
function and variational parameters aq and j3q all real. 
Notice that aq controls the extent of the Gaussian wave 
function, i.e., the amplitude of lattice vibration, and that 


1 {Qq — Pq)'^ 
^ 2 al 
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FIG. 2. (color online) Ground state phase diagram of the 
two-band Hubbard model in (a) U-U' plane {gji = 0.19) and 
(b) U-g plane {U' = P/2). We set Ue = Uf = U and uj/t = 
0.1. NM, CDWI, and BI denote normal metal, charge-density- 
wave insulator, and band insulator, respectively. The electron 
density is n = 2, i.e., at half filling. 


Pq corresponds to the average value of Qq and thus there 
exists a static lattice distortion with {up 0 for a finite 
Pq. The Monte Carlo update scheme for {Qq} and the 
estimation of phonon energy are the same as in Ref. |4^ . 
The remaining part is an electron-phonon projection 


operator: 


Pe-ph = exp 





This oper¬ 


ator controls the attraction between c electrons and / 
holes which results from the electron-phonon interaction 
and 7 is a variational parameter. 

The variational parameters in jT) are therefore tc, t}, 
He, Vq, a, {gr}, [aq], [Pq], and 7 , and they 

are simultaneously optimized using stochastic reconfig¬ 
uration method |43l |. The system sizes are varied from 
L X L = 12 X 12 to 24 X 24 with antiperiodic boundary 
conditions in both directions of primitive lattice vectors 
for the triangular lattice. 


III. RESULT 

Figure (Ha) shows the ground state phase diagram 
where Ue = Uf = U and U' are varied for fixed 
g{k,q)/t = g/t = 0.19 and u}{q)/t = uj/t = 0.1 [i^j . 
We find that there are three distinct phases in the phase 
diagram: normal metal (NM), charge-density-wave in¬ 
sulator (CDWI), and band insulator (BI). When U' is 
large enough, the c band is lifted above the Fermi energy 
and the BI phase with the empty c band and the fully- 
occupied / band is stabilized. No static lattice distortion 
is observed in both NM and BI phases. 

Between the NM and BI phases, the CDWI phase 
emerges where the c-f hybridization parameter 

” ^^i^k+qafko- + H.C.) (7) 

k.<7 

is finite for q corresponding to the three ordering 
wave vectors qi, q 2 , and qg, simultaneously, implying a 
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FIG. 3. (color online) (a) Average c-electron density (rix) (red 
bars) and /-hole density 2 — (n^) (blue bars) in 2 x 2 unit 
cell for X = A, B, C, and D [see Fig. [T][c)]. (b) Snap shots of 
c-f bond length ua and mb at A and B sites, respectively, as 
a function of Monte Carlo (MC) step. The model parameters 
used are {U/t,U'/t, g/t,u}/t) = {4:.0, 2.0, 0.19, 0.1) for L = 24 
in the CDWI phase. 


FIG. 4. (color online) Momentum resolved c-f hybridization 
4>{k) for (a) {U/t,U'/t)=i8.0, 4.0) and (b) iU/t,U'/t)={3.0, 
1.5). The noninteracting Fermi momentum fep (folded around 
F point at the center) and fep are shown with white solid and 
black dashed curves, respectively [see also Fig. [T][a)[. gft = 
0.19 and uj/t = 0.1 are fixed for L — 24. 


triple-q CDW state. Here, (O) = (\E'|C>|'I')/(4'|tE') with 
the optimized |T'). Thus, the first Brillouin zone is folded 
as indicated in Fig. (Ha). It leads to the charge dispro¬ 
portionation in 2x2 unit cell [Fig. [TJc)] with one charge 
rich A site and three charge poor B, C, and D sites [see 
Fig. [31[a)]. In the CDWI phase, a static lattice distor¬ 
tion always occurs through the electron-phonon interac¬ 
tion and the “pure” exciton condensation without lattice 
distortion is never found. Note that the CDWI phase 
is limited to a narrow region in Fig. [H[a) especially for 
small U/t in spite of a finite g/t. This is in sharp contrast 
with the case of a square lattice where the NM phase ap¬ 
pears only at [/' = 0 and the exciton condensation phase 
is widely observed even without the electron-phonon in¬ 
teraction [H, 113, 13 . This difference is caused by the 
different FS nesting condition: the FS nesting is better 
(perfect if only with the nearest-neighbor hopping) in the 
square lattice but poor in the triangular lattice. There¬ 
fore, the electron-phonon interaction is indispensable to 
manifest the CDWI phase under the poor FS nesting 0 . 

We also show the phase diagram in Fig. [2][b) where 
U and g are varied with U' = U/2. The CDWI re¬ 
gion is enlarged with increasing U and g, implying that 
both Coulomb interaction and electron-phonon interac¬ 
tion stabilize the CDWI phase. This result is thus qual¬ 
itatively consistent with previous study [l^. We also 
find that the CDWI phase is not stabilized, but only 
the NM and BI phases appear, when g = 0, at least, 
in a realistic parameter region. This suggests that the 
“pure” exciton condensation induced by the Coulomb 
interaction alone, the original idea of exciton conden¬ 
sation [2§ - [30j| . is difficult to realize in our model. The 
pure exciton condensation certainly occurs in particular 
models such as one-dimensional models or two- 

dimensional models with perfectly nested electron and 
hole FSs m, [13, m, [53 ■ Therefore, the stability of the 


pure exciton condensation depends strongly on the lat¬ 
tice structure and the underlying FS. 

Let us now examine the detailed properties of the 
CDWI phase. Figure [3][a) shows the distribution of aver¬ 
age c-electron density (n^) and /-hole density 2 — (n^) in 
2x2 unit cell (X=A, B, C, and D). It is found in Fig.[^a) 
that c electrons and / holes, i.e., mobile carriers, are con¬ 
centrated mostly at A site with {n'^) (n^) > 2 to gain 

the c-f hybridization energy coupled with the lattice dis¬ 
tortion, while the number of these mobile carriers are 
small and (n^) + (n^) < 2 at B, C, and D sites. There¬ 
fore, the system clearly exhibits the charge dispropor¬ 
tionation. Notice that the mobile carrier densities at B, 
C, and D sites are the same within the statistical errors 
simply because the three ordering wave vectors qi, q 2 , 
and 93 are equivalent in a hexagonal lattice structure. 

Next, let us discuss the lattice degrees of freedom in the 
CDWI phase. In the VMC calculation, the bond length 
Ui always fluctuates around the average value during the 
Monte Carlo steps [IJ]. As shown in Fig. EKb), we find 
that the bond length at A site is shortened from the 
original one (ma < 0), while those at B, C, and D sites 
are lengthened (only ub > 0 is shown). Figs. EJa) and 
[31[b) thus confirm that the electronic and lattice degrees 
of freedom are strongly coupled and cooperatively induce 
the CDWI phase. 

Furthermore, we study the character of electron-hole 
pairing from the viewpoint of BCS-BEC crossover, which 
has been often discussed in the exciton problems [^, 148 . 
[sol - fs^ . For this purpose, we calculate the momentum 
resolved c-f hybridization (/)(fc) defined as 

</>(fe)=E(4+,.^<^ + H.c.). (8) 

q,a 

Figure m shows (()(fc) for {U/t,U'/t)={8.0, 4.0) and (3.0, 
1.5), both being located in the CDWI phase in Fig. [21 
For {U/t, [/'/t)=(8.0, 4.0), 4>{k) is extended in the whole 
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Brillouin zone, indicating the strong-coupling BEC-like 
pairing due to the large Coulomb interaction. Although 
4>{k) becomes less extended with decreasing the Coulomb 
interaction, it still has a broad structure away from 
the Fermi momentum fcp, as shown in Fig. |4l[b) for 
{U/t,U'/t)={3.0, 1.5). Indeed, the CDWI region rapidly 
decreases with decreasing U [see Fig. [2][a)] and our sys¬ 
tematic calculations do not find a clear BCS-like region in 
the CDWI phase shown in Fig. [2l Because of i) the poor 
nesting between c-electron and /-hole FSs and ii) the 
small density of states around the Fermi energy for low 
carrier densities, the energy gain due to the gap open¬ 
ing induced by the c-/ hybridization in the vicinity of 
fcp is small and hence the weak-coupling BCS-like pair¬ 
ing is not favored. Even in such a case, the BEC-like 
tightly-bounded electron-hole pairing in real space can 
be induced by the electron-phonon interaction with the 
help of Coulomb interaction and dominates the CDWI 
phase. In contrast, we have found a clear and wide BCS- 
like region for the same model but in a square lattice 
with perfectly nested FSs . Therefore, the FS nesting 
is essential for the BCS-like pairing. 


IV. DISCUSSION 

Finally, let us discuss the implication of our results 
for IT-TiSe 2 . In our model, Ti 3d and Se 4p bands are 
simplified as c and / bands, respectively, and the or¬ 
bital characters are ignored. Moreover, our model only 
includes the change of c-f bond length which couples 
with exciton condensation. Even with these simplihca- 
tions, our model captures the important energy scales 
of IT-TiSe 2 . The electron-phonon coupling used here is 
ju) = 0 — 3.6t Ri 0 — 1.8 eV, which is relevant for 
lT-TiSe 2 [ll if we take t r; 0.5 eV. The lattice distor¬ 
tion obtained in our calculation is 0.05 — 0.2A, consistent 
with the observed value ~ 0.085A sm. Our results 
thus suggest that the CDW phase observed in ir-TiSe 2 
is due to the strong-coupling BEC-like electron-hole pair¬ 
ing. Indeed, the BEC-like character is indicated by sev¬ 
eral theoretical works [mis and experimental obser¬ 
vations such as a short coherence length estimated by 
Kohn anomaly @ , lack of incommensurate CDW phase, 
relatively high electrical resistivity above Tc, and a large 
value of 2A/fcBT'c (A: the CDW gap) [s^. 

On the other hand, the chiral CDW phase observed 
in IT-TiSe 2 [H, [5^ is beyond our model. In the chi¬ 
ral CDW phase, the charge density is modulated with 
clockwise or anticlockwise pattern. The proper descrip¬ 


tion of this phase requires three dimensionality or 
higher-order electron-phonon and phonon-phonon inter¬ 
actions [ 2 ^ which induce the phase difference between 
the three ordering wave vectors. The origin of the SC in¬ 
duced by applying pressure or intercalation of Cu atoms 
is also an interesting unresolved issue. The relation be¬ 
tween the CDW and the SC is still controversial [H, 
and both conventional [ 2 O, [ 2 ^ and unconventional 
SC [m have been proposed. Our results suggest that 
both electronic and lattice degrees of freedom are crucial 
to understand the origin of the SC. Our study will be 
a first step toward the unified understanding of various 
quantum phases observed in IT-TiSe 2 . 

V. SUMMARY 

In summary, we have studied the two-band Hubbard 
model in a triangular lattice for IT-TiSe 2 with the 
electron-phonon interaction. The VMC method is em¬ 
ployed to treat the electronic and lattice degrees of free¬ 
dom on an equal footing beyond the mean-field approxi¬ 
mation. We have shown that both Coulomb and electron- 
phonon interactions stabilize the CDW phase. We have 
found that the “pure” exciton condensation without the 
lattice distortion is difficult to realize and the electron- 
phonon interaction is essential for the CDW phase. The 
character of electron-hole pairing within the CDW phase 
has also been examined by calculating the momentum re¬ 
solved c-f hybridization. We have shown that the strong¬ 
coupling BEC-like pairing dominates the CDW phase. 
Under the poor FS nesting condition and with small den¬ 
sity of states around Fermi energy, the energy gain due 
to the gap opening in the vicinity of fep is small and 
hence the weak-coupling BCS-like pairing is not favored. 
Our results thus conclude that the CDW phase observed 
in IT-TiSe 2 originates from the strong-coupling BEC-like 
electron-hole pairing due to the cooperative Coulomb and 
electron-phonon interactions. 

ACKNOWLEDGMENT 

The authors thank Y. Fuseya, T. Shirakawa, T. 
Kaneko, and K. Imura for useful discussions. The compu¬ 
tation has been done using the RIKEN Cluster of Clus¬ 
ters (RICC) facility and the facilities of the Supercom¬ 
puter Center, Institute for Solid State Physics, Univer¬ 
sity of Tokyo. This work has been supported by JSPS 
KAKENHI Grant No. 26800198 and in part by RIKEN 
iTHES Project. 


[1] G. Griiner, Density Waves in Solids (Perseus Publishing, 
New York, 2000). 


[2] M. D. Johannes and I. I. Mazin, Phys. Rev. B 77, 165135 
(2008). 



6 


[3] J. A. Wilson, F. J. Di Salvo, and S. Mahajan, Phys. Rev. 
Lett. 32, 882 (1974). 

[4] F. J. Di Salvo, D. E. Moncton, and J. V. Waszczak, Phys. 
Rev. B 14, 4321 (1976). 

[5] K. Rossnagel, J. Phys.: Condens. Matter 23, 213001 

( 2011 ). 

[6] B. Sipos, A. F. Kusmartseva, A. Akrap, H. Berger, L. 
Forro, and E. Tutis, Nat. Mater. 7, 960 (2008). 

[7] A. F. Kusmartseva, B. Sipos, H. Berger, L. Forro, and E. 
Tutis, Phys. Rev. Lett. 103, 236401 (2009). 

[8] L. J. Li, W. J. Lu, X. D. Zhu, L. S. Ling, Z. Qn, and Y. 
P. Sun, Europhys. Lett. 97, 67005 (2012). 

[9] Y. Liu, R. Ang, W. J. Lu, W. H. Song, L. J. Li, and Y. 
P. Sun, Appl. Phys. Lett. 102, 192602 (2013). 

[10] E. Morosan, H. W. Zandbergen, B. S. Dennis, J. W. G. 
Bos, Y. Onose, T. Klimczuk, A. P. Ramirez, N. P. Ong 
and R. J. Cava, Nat. Phys. 2, 544 (2006). 

[11] G. Li, W. Z. Hu, D. Qian, D. Hsieh, M. Z. Hasan, E. 
Morosan, R. J. Cava, and N. L. Wang, Phys. Rev. Lett. 
99, 027404 (2007). 

[12] J. C. E. Rasch, T. Stemmier, B. Muller, L. Dudy, and R. 
Manzke, Phys. Rev. Lett. 101, 237602 (2008). 

[13] G. Monney, C. Monney, B. Hildebrand, P. Aebi, and H. 
Beck, Phys. Rev. Lett. 114, 086402 (2015). 

[14] A. Zunger and A. J. Freeman, Phys. Rev. B 17, 1839 
(1978). 

[15] J. A. Wilson, Solid State Commun. 22, 551 (1977). 

[16] H. Cercellier, C. Monney, F. Clerc, C. Battaglia, L. De- 
spont, M. G. Gamier, H. Beck, P. Aebi, L. Patthey, H. 
Berger, and L. Forro, Phys. Rev. Lett. 99, 146403 (2007). 

[17] M. Cazzaniga, H. Cercellier, M. Holzmann, C. Monney, 
P. Aebi, G. Onida, and V. Olevano Phys. Rev. B 85, 
195111 (2012). 

[18] H. P. Hughes, J. Phys. C: Solid State Phys. 10, L319 
(1977). 

[19] N. Suzuki, A. Yamamoto, and K. Motizuki, J. Phys. Soc. 
Jpn. 54, 4668 (1985). 

[20] M. Calandra and F. Mauri, Phys. Rev. Lett. 106, 196406 

( 2011 ). 

[21] Z. Zhu, Y. Cheng, and U. Schwingenschlogl, Phys. Rev. 
B 85, 245133 (2012). 

[22] F. Weber, S. Rosenkranz, J.-P. Castellan, R. Osborn, 
G. Karapetrov, R. Hott, R. Held, K.-P. Bohnen, and A. 
Alatas, Phys. Rev. Lett. 107, 266401 (2011). 

[23] J. van Wezel, P. Nahai-Williamson, and S. S. Saxena, 
Phys. Rev. B 83, 024502 (2011). 

[24] A. Taraphder, S. Koley, N. S. Vidhyadhiraja, and M. S. 
Laad, Phys. Rev. Lett. 106, 236405 (2011). 

[25] C. Monney, G. Monney, P. Aebi, and H. Beck, New J. 
Phys. 14, 075026 (2012). 

[26] B. Zenker, H. Fehske, H. Beck, C. Monney, and A. R. 
Bishop, Phys. Rev. B 88, 075138 (2013). 

[27] S. Koley, M. S. Laad, N. S. Vidhyadhiraja, and A. Tara¬ 
phder, Phys. Rev. B 90, 115146 (2014). 

[28] N. F. Mott, Phil. Mag. 6 , 287 (1961). 

[29] A. N. Kozlov and L. A. Maksimov, Sov. Phys. JETP 21, 
790 (1965). 

[30] D. Jerome, T. M. Rice, and W. Kohn, Phys. Rev. 158, 
462 (1967). 

[31] B. 1. Halperin and T. M. Rice, Rev. Mod. Phys. 40, 755 
(1968). 

[32] B. Bucher, P. Steiner, and P. Wachter, Phys. Rev. Lett. 
67, 2717 (1991). 


[33] T. Mizuno, T. lizuka, S. Kimura, K. Matsubayashi, K. 
Imura, H. S. Suzuki, and N. K. Sato, J. Phys. Soc. Jpn. 
77, 113704 (2008). 

[34] Y. Wakisaka, T. Sudayama, K. Takubo, T. Mizokawa, 
M. Arita, H. Namatame, M. Taniguchi, N. Katayama, 
M. Nohara, and H. Takagi, Phys. Rev. Lett. 103, 026402 
(2009). 

[35] K. Seki, Y. Wakisaka, T. Kaneko, T. Toriyama, T. Kon- 
ishi, T. Sudayama, N. L. Saini, M. Arita, H. Namatame, 
M. Taniguchi, N. Katayama, M. Nohara, H. Takagi, T. 
Mizokawa, and Y. Ohta, Phys. Rev. B 90, 155116 (2014). 

[36] J. Zittartz, Phys. Rev. 162, 752 (1967). 

[37] H. Watanabe, T. Shirakawa, and S. Yunoki, Phys. Rev. 
Lett. 105, 216410 (2010); H. Watanabe, T. Shirakawa, 
and S. Yunoki, Phys. Rev. Lett. 110, 027002 (2013); H. 
Watanabe, T. Shirakawa, and S. Yunoki, Phys. Rev. B 
89, 165115 (2014). 

[38] C. Monney, E. F. Schwier, M. G. Gamier, N. Mariotti, C. 
Didiot, H. Cercellier, J. Marcus, H. Berger, A. N. Titov, 
H. Beck, and P. Aebi, New J. Phys. 12, 125019 (2010). 

[39] H. Watanabe, K. Seki, and S. Yunoki, J. Phys.: Conf. 
Ser. 592, 012097 (2015). 

[40] J. Biinemann, W. Weber, and F. Gebhard, Phys. Rev. B 
57, 6896 (1998). 

[41] B. J. Alder, K. J. Runge, and R. T. Scalettar, Phys. Rev. 
Lett. 79, 3022 (1997). 

[42] T. Ohgoe and M. Imada, Phys. Rev. B 89, 195139 (2014). 

[43] S. Sorella, Phys. Rev. B 64, 024512 (2001); S. Yunoki 
and S. Sorella, Phys. Rev. B 74, 014408 (2006). 

[44] The system size used is 16 x 16. However, we have checked 
that the phase boundary is almost unchanged for L > 16. 

[45] Aq is the order parameter of exciton condensation in the 
absence of lattice distortion. However, as pointed out by 
Zenker et al. [i^, finite Aq does not necessarily corre¬ 
spond to the exciton “condensation” in the presence of 
lattice distortion due to electron-phonon interactions, in 
which the phase coherence of exciton pairs, for example, 
should be discussed. 

[46] B. Zenker, H. Fehske, and H. Beck, Phys. Rev. B 90, 
195118 (2014). 

[47] B. Zocher, C.Timm, and P. M. R. Brydon, Phys. Rev. B 
84, 144425 (2011). 

[48] T. Kaneko, K. Seki, and Y. Ohta, Phys. Rev. B 85, 
165135 (2012). 

[49] T. Kaneko, T. Toriyama, T. Konishi, and Y. Ohta, Phys. 
Rev. B 87, 035121 (2013). 

[50] K. Seki, R. Eder, and Y. Ohta, Phys. Rev. B 84, 245106 

( 2011 ). 

[51] P. Nozieres and S. Schmitt-Rink, J. Low. Temp. Phys. 
59, 195 (1985). 

[52] F. X. Bronold and H. Fehske, Phys. Rev. B 74, 165107 
(2006). 

[53] B. Zenker, D. Ihle, F. X. Bronold, and H. Fehske, Phys. 
Rev. B 85, 121102(R) (2012). 

[54] In our model, only two bands related to the exciton con¬ 
densation are considered and the resulting CDW phase 
is an insulator. In ir-TiSe 2 , the system remains metal¬ 
lic even when the CDW transition occurs. This is due to 
the imperfect opening of the gap. Therefore, the lattice 
distortion obtained in our model tends to be larger than 
the experimental one. 

[55] J. Ishioka, Y. H. Liu, K. Shimatake, T. Kurosawa, K. 
Ichimura, Y. Toda, M. Oda, and S. Tanda, Phys. Rev. 
Lett. 105, 176401 (2010). 


7 


[56] J.-P. Castellan, S. Rosenkranz, R. Osborn, Q. Li, K. E. 
Gray, X. Luo, U. Welp, G. Karapetrov, J. P. C. Ruff, and 
J. van Wezel, Phys. Rev. Lett. 110, 196404 (2013). 

[57] J. van Wezel, Europhys. Lett. 96, 67011 (2011). 

[58] D. Qian, D. Hsieh, L. Wray, E. Morosan, N. L. Wang, Y. 
Xia, R. J. Cava, and M. Z. Hasan, Phys. Rev. Lett. 98, 
117007 (2007). 


[59] Y. 1. Joe, X. M. Chen, P. Ghaemi, K. D. Finkelstein, G. 
A. de la Pena, Y. Gan, J. C. T. Lee, S. Yuan, J. Geek, G. 
J. MacDougall, T. C. Chiang, S. L. Cooper, E. Fradkin 
and P. Abbamonte, Nature Phys. 10, 421 (2014). 

[60] S. Y. Li, G. Wu, X. H. Chen, and L. Taillefer, Phys. Rev. 
Lett. 99, 107001 (2007). 

[61] R. Ganesh, G. Baskaran, J. van den Brink, and D. V. 
Efremov, Phys. Rev. Lett. 113, 177001 (2014). 


